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The noncommutative spectral action extends our familiar notion of commutative spaces, using the 
data encoded in a spectral triple on an almost commutative space. Varying a rather simple action, 
one can derive all of the standard model of particle physics in this setting, in addition to a modified 
version of Einstein-Hilbert gravity. In this letter we use observations of pulsar timings, assuming 
that no deviation from General Relativity has been observed, to constrain the gravitational sector of 
this theory. Whilst the bounds on the coupling constants remain rather weak, they are comparable 
to existing bounds on deviations from General Relativity in other settings and are likely to be further 
constrained by future observations. 

PACS numbers: ll.10.Nx, 04.50.+h, 12.10.-g, 11.15.-q, 12.10.Dm 



INTRODUCTION 



Approaching Planckian energies, the assumption of 
Riemannian geometry and the validity of General Rel- 
ativity (GR) break down and one is forced to describe 
the space-time geometry within a fully quantum frame- 
work. NonCommutative Geometry (NCG) [l|, |2| is based 
on the idea that as we approach Planckian energy scales, 
our intuitive description of space-time being a commuta- 
tive manifold ceases to be a valid approximation. In its 
simplest but nevertheless powerful version, NCG implies 
that slightly below Planck energy, space-time is well ap- 
proximated by the product of a four-dimensional smooth 
compact Riemannian manifold A4 and a finite noncom- 
mutative space J ' . Such spaces are called "almost com- 
mutative" spaces and they are the simplest extensions of 
the commutative spaces we use in GR. This is a strong 
assumption which is certainly expected to break at the 
Planck scale, where the notion of classical geometry loses 
all meaning, however at low energies it should be a good 
approximation. 

The noncommutative nature of J- is given by the real 
spectral triple (A,ri,D), where A is an involution of op- 
erators on the Hilbert space H, and D is a self-adjoint 
unbounded operator in W. The algebra A is the algebra 
of coordinates, the operator D corresponds to the inverse 
line element of Riemannian geometry, and the commuta- 
tor [D, a] with a E A plays the role of the differential 
quotient da/ds, with ds the unit of length. The result- 
ing physical Lagrangian is obtained from the asymptotic 
expansion in the energy scale A of a spectral action func- 
tional of the form Tr(/(D/A)) defined on noncommuta- 
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tive spaces, where / is a cut-off function (i.e. a test func- 
tion of compact support). The coupling with fermions 
can be obtained by including an additional term in the 
spectral action functional. The choice of the finite dimen- 
sional algebra is the underlying geometric input which 
determines the physical implications of the model, in par- 
ticular the particle content of the theory. 

The NCG spectral action offers a simple and ele- 
gant explanation for the phenomenology of the Stan- 
dard Model (SM) compatible with right-handed neutri- 
nos and neutrino masses |3[ (the approach has also been 
used to derive supersymmetric extensions to the standard 
model [4J]). This approach to the SM has been proposed 
as a way to achieve unification, based on the symplectic- 
unitary group (the algebra constructed in M. x T is as- 
sumed to be symplectic-unitary) in the Hilbert space, in- 
stead of the finite dimensional Lie groups. Note that the 
NCG spectral action is a classical theory which, in prin- 
ciple, should eventually be quantized. Whilst an under- 
standing of how to quantize such noncommutative spaces 
has not yet been fully developed, already at the classical 
level the theory introduces several extensions to standard 
GR. Specifically, the physical Lagrangian contains, in ad- 
dition to the full SM Lagrangian, the Einstein-Hilbert ac- 
tion with a cosmological term, a topological term related 
to the Euler characteristic of the space-time manifold, 
a conformal Weyl term and a conformal coupling of the 
Higgs field to gravity. In contrast to the SM on a fixed 
background, the coefficients of the gravitational terms in 
this NCG action depend on the Yukawa parameters of 
the particle physics content. 

The parameters of the NCG spectral action model are 
set at the scale A, considered to be the unification scale, 
while physical predictions at lower energies are recovered 
by running the parameters down through Renormaliza- 
tion Group Equations (RGE). Thus, the spectral action 
at the unification scale A is directly applicable to early 
universe cosmological models [a-ISl, while extrapolations 



to lower energies can be obtained via RGE and inclusion 
of nonperturbative effects in the spectral action. 

The gravitational part of the asymptotic formula for 
the bosonic sector of the NCG spectral action, includ- 
ing the coupling between the Higgs field and the Ricci 
curvature scalar, reads |3| 
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H is a rescaling H = (%/ a /o/ 7r )0 of the Higgs held <f> 
to normalize the kinetic energy, the momentum /o is 
physically related to the coupling constants at unifica- 
tion and the coefficient a is related to the fermion and 
lepton masses and lepton mixing. Note that we are us- 
ing conventions in which the signature is (—,+,+,+) 
and the Ricci tensor is defined as R^ v = R p p, up , with 
R^ P a ^a = [\/V> Vf] w p- In the above action, Eq. (JTJ, 
the first term is the familiar Einstein-Hilbert term, the 
second one is a Weyl curvature term, the third term 
R*R* = (1/ '4)e^P a 'e aPlS R^Rf a , is the topological term 
that integrates to the Euler characteristic, hence nondy- 
namical, and the last one couples gravity to the Higgs 
field 1 and can have important consequences at high en- 
ergies, such as in the early universe [5-8]. Here, we will 
be concerned with the low energy, weak curvature regime 
where this term is small. 

Neglecting the nonminimal coupling between the Higgs 
field and the Ricci curvature, the equations of motion 
derived from the spectral action above read [5j 
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where (3 2 is defined as (3 2 = — I /(327rGao). Notice in par- 
ticular that the NCG corrections vanish for Fricdmann- 
Lemaitre- Robertson- Walker (FLRW) cosmologies [5| and 
Schwarzschild solutions, which makes it difficult to place 
restrictions on these terms via cosmology or solar-system 
tests. The best constraint on, different ad hoc, curva- 
ture squared terms are obtained from measurements of 
the orbital precession of Mercury, imposing a rather weak 
lower bound on (3, namely (3 > 3.2 x I0~ 9 m _1 [9]. This 
constraint was however found for terms of different form 
(but of the same order) to the Weyl term appearing in 
the NCG spectral action approach we investigate here. 
In what follows, we will specifically study how one can 
constrain (3 within the NCG context. The parameter /3 
can be equivalently expressed in terms of /o, through 



1 Such a term should always be present when one considers gravity 
coupled to scalar fields. 



(3 2 = (57r)/(48G/o), so by imposing a lower limit to (3, 
we actually set an upper limit to the moment /o of the 
cut-off function used to define the spectral action. The 
normalization of kinetic terms in the spectral action im- 
poses the following relation between the gauge couplings 
of the Standard Model, gi,g2, 33 and the coefficient /o [3(] , 
namely <? 3 2 /o/(2^ 2 ) = 1/4 , g 2 = g 2 = (5/3) 5 2 . The im- 
portance of constraining f3 is thus clear, since /o can be 
used to specify the initial conditions on the gauge cou- 
plings, a constraint on (3 corresponds to a restriction on 
the particle physics at unification. 

We will study the energy lost to gravitational radia- 
tion by orbiting binaries, so we consider the weak field 
limit of Eq. ([2]). The general first order solution for a 
perturbation against a Minkowski background is [I0( 
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T = t— if is the difference between the time of observation 
(t) and emission (£') of the perturbation, R = r — r' is the 
difference between the locations of the observer (r) and 
emitter (r'), J\ is a Bessel function of the first kind and 
is the Heavyside step function. In the far field limit, 
|r| sa |r — r'|, the spatial components of Eq. ^ become 
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where we have, introduced the quadrupole moment, 

D lk (t) = ^ / dr a;^ fe r 00 (r, t) . (5) 

From Eq. ([2]) is it clear that this theory reduces to that 
of GR in the (3 — > oo limit, and one can check that in 
this limit Eq. ^ does indeed reproduce the standard 
result for a massless graviton. For finite f3 however, one 
finds that gravitational radiation contains both massive 
and massless modes, both of which are sourced from the 
quadrupole moment of the system. 



GRAVITATIONAL RADIATION FROM 
CIRCULAR BINARIES 

We will derive the explicit formula for the energy lost 
to gravitational radiation from a binary pair in a circu- 
lar orbit. One can similarly consider binaries in elliptical 
orbits; for simplicity we consider only circular Keplerian 
orbits. Similarly, we neglect effects due to the internal 



structure of the bodies as well as local astrophysical ef- 
fects, such as mass transfer, tidal stripping etc., focusing 
instead on the purely gravitational consequences of NCG. 
Consider a circular binary pair, of masses mi, mi. For 
such a system, orbiting in the xy-plane, the only nonzero 
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components of the quadrupole moment are [11 



D xx (t) = 12p\p\ 2 sin (2^ (t))uj 3 

= -D yy (t) , 
D xv (t) = -I2p\p\ 2 cos(2tp(t))oj 3 , 
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where /i = m\raij(rn\ + 7712) is the reduced mass of the 
system, \p\ is the magnitude of the separation vector be- 
tween the bodies, which is constant for circular orbits, ip 
is the angle of the bodies relative to the x-axis and lj = %jj 
is the orbital frequency, which for this simple system is a 
constant given by 



u = ip = |p| -3/ V<3(mi + m 2 ) 
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Following the standard approach (see e.g., Ref. [Hi), the 
rate of energy loss , in the far field limit, is 
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This allows us to explicitly test the theory by compar- 
ing this prediction to binary pulsar measurements, for 
which the energy loss has been very well characterized 
(see, Table H| and hence constrain /3. 
Using Eq. (0} one finds [HI 
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where we have defined the functions 
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The integrals above, Eq. (TIU)) . exhibit a strong resonance 
behavior at z = 1, however they are easily evaluated for 
both z < 1 and z > 1. This resonance corresponds to a 
critical frequency given by 
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and we can expect strong deviations from the standard 
results of GR for orbital frequencies close to this critical 
frequency. 

One can evaluate numerically the functions in Eq. (jTDJ) 
and fit them to an explicit functional form. Thus, for 
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where C is approximately a constant, C w 0.175, except 
as oj approaches uj c . Similarly, for uj > uj c one gets 
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where the function / is approximately 
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Note that as we will see, the precise form of this function 
is unimportant. See Ref. [lfj for a discussion on the 
accuracy of the approximations. 

Using the above approximations, one can expand 
Eq. © in the large distance (large |r|) limit, to find the 
rate of energy lost to gravitational radiation: 

32Gp 2 p 4 uj 6 




5c 5 
c 



;Ji /3M 



4sin 2 (/3|r|/(^)) 



;w < u} c 



,(16) 



where in the u> < u c case the dots refer to higher pow- 
ers of 1/ (/3|r|). Thus, for orbital frequencies small com- 
pared to u> c , any deviation from the standard result is 
suppressed by the distance to the source. Notice that in 
this case, the f3 — > 00 (i.e., ao — > 0) limit reproduces the 
GR result, as it should. For the u) > u> c case, the result 
would only agree with the General Relativistic result if 
j3\r\f (u/ui c ) — tt/3, which is clearly not true for systems 
at different distances, |r|, with different orbital frequen- 
cies, ui. Hence, we can immediately eliminate the u > u c 
case, simply by noting that observations of the energy 
lost to gravitational radiation agree, to a high level of 
accuracy, with those of GR for many different systems. 

The resonance appearing in Eq. ([TO]) leads to a sim- 
ple heuristic argument to rule out the ui > u c case. A 
system with uj < cj c cannot increase its orbital frequency 
above w c , without losing a significant fraction of its en- 
ergy to gravitational radiation. Similarly, a system with 
uj > uj c cannot decrease its orbital frequency across this 



boundary. Since one expects all astrophysical systems 
to have formed from the coalescence of relatively cold, 
slowly moving systems, it is reasonable to suppose that at 
some time in the past, all binary systems had very slowly 
varying quadrupole moments and hence that lu < u) c . In 
the following, we will see that this restriction places a 
strong constraint on (which defines oj c ). 

For the physically interesting case of u < lu C7 the am- 
plitude of the deviation from the standard result is small, 
due to the l/|r| suppression, however there are two in- 
teresting features: firstly, the existence of a critical fre- 
quency uj c and secondly, the fact that the rate of flux of 
gravitational radiation is oscillatory. 

The critical frequency comes from the fact that this 
theory contains a natural frequency scale given by 0c ~ 
c (— aoG)~ . This scale is set by the moments of the cut- 
off function used to define the spectral action. Physically 
one can think of this as the scale at which noncommuta- 
tivc effects become dominant. What is important for this 
work, is that the binary systems must have orbital fre- 
quencies below this critical value, since otherwise the the- 
ory would predict significant deviations from the results 
of standard GR, which can be ruled out observationally. 

The presence of the Bessel function in Eq. (|16|) means 
that the amplitude of the deviation from the standard 
result of GR will oscillate both with changing distances 
and changing frequencies. Whilst such correlations may 
present new observational signatures, the effect is heavily 
suppressed by the |r | 1 factor in Eq. (fT6|) : here we focus 
on overall amplitude of deviations from the GR result. 



ASTROPHYSICAL CONSTRAINTS 

Having calculated the general form of the gravitational 
radiation from binary systems, within this NCG theory 
of gravity, we can now constrain the main parameter of 
the theory, 0, via observational data. From Eq. (fT6| it is 
clear that the only data needed is the orbital frequency 
and the distance to the binary system. If one had consid- 
ered elliptical binary orbits, additional parameters would 
come into play, however we are concerned only with the 
order of magnitude of the constraint and hence neglect 
such additional complications. A more detailed quantita- 
tive analysis would require the inclusion of the ellipticity 
as well as various other near field effects. 

Table U gives the binaries we are considering. We focus 
on binary pulsars for which the rate of change of the or- 
bital frequency has been well characterized. In all these 
cases the predictions of GR agree with the data to high 
accuracy (see, Table QJ. We can thus restrict by requir- 
ing that the magnitude of deviations from GR, given by 
Eq. (fl"6|) . be less than this uncertainty. 

Using the data on the six binaries given in Table HI and 
requiring that uj < lu c (see, the discussion above), we find 
> 7.55 x 10~ 13 m^ 1 . The restrictions coming from the 



Binary 


Distance 


Orbital 


Eccentricity 


GR 




(PC) 


Period (hr) 




(%) 


PSR J0737-3039 


- 500 


2.454 


0.088 


0.2 


PSR JI0I2-5307 


- 840 


14.5 


< 10~ 6 


10 


PSR JII4I-6545 


> 3700 


4.74 


0.17 


6 


PSRBI9I6+I6 


- 6400 


7.752 


0.617 


0.1 


PSR BI534+I2 


- 1100 


10.1 


? 


1 


PSR B2127+11C 


- 9980 


8.045 


0.68 


3 



TABLE I: We calculate the constraint on the NCG 
theory of gravity, via the predicted energy lost 
to gravitational radiation from the above binaries 
(Refs. 0, IH, 0, IH, 0. E3l> respectively). The 
column marked GR, indicates the approximate accuracy to 
which the rate of change of the orbital period agrees with 
the predictions of GR. 



PSR J0737-3039 


> 7.55 x 10~ 13 m- 1 


PSR J1012-5307 


> 7.94 x 10~ 14 m- 1 


PSR J1141-6545 


> 3.90 x 10~ 13 m- 1 


PSR B1913+16 


> 2.39 x 10~ 13 m- 1 


PSR B1534+12 


> 1.83 x 10~ 13 iri -1 


PSR B2127+11C 


> 2.30 x 10~ 13 m" 1 



TABLE II: For each binary system, we restrict by requiring 
that the energy lost to gravitational radiation agrees with the 
prediction of GR to within observational uncertainties. 



individual systems are given in Table |TTJ 

Due to the large distances to these systems, the con- 
straint is almost exactly due to w < u c which, using the 
definition of uj c given in Eq. (fTTj). becomes > 2u/c. 
Thus, the strongest constraint comes from systems with 
high orbital frequencies. This will be true for all systems 
for which 2w|r|/c is large. Future observations of rapidly 
orbiting binaries, relatively close to the Earth, could thus 
improve this constraint by many orders of magnitude. 

This dependence of the constraint on the orbital fre- 
quency, suggests that other astrophysical objects, with 
high frequency periodicity, such as individual pulsars 
or merger in-spirals may provide a more stringent con- 
straint. Whilst the analysis given here is only applicable 
to binaries, it can be extended by replacing Eq. © by the 
quadrupole moments of whatever system is of interest. 



CONCLUSIONS 

General Relativity is formulated within the arena of 
Riemannian geometry a natural extension of which is 
NonCommutative Geometry. The spectral action ap- 
proach produces all the Standard Model fields as well 
as gravitational terms, from purely geometric consider- 
ations. Thus, both gravity and matter are treated in a 



similar manner within NCG, which also provides us with 
concrete relationships between matter and gravitational 
couplings. The asymptotic expansion of the gravitational 
sector of this theory produced modifications to GR and 
in this paper we use these modifications to test and con- 
strain the theory through observations. 

We have considered the energy lost by circular binary 
systems to gravitational radiation and shown that for 
the predicted values to agree with observations, a key 
parameter of the theory can be constrained. We have 
focused on binary pulsar systems, for which the rate of 
change of the orbital frequency is well known and ex- 
plicitly calculated the predicted deviation from the GR 
expressions. We have shown that this restricts the value 
of the Weyl squared coupling in the bosonic action (i.e., 
ao hi Eq. ([!])). This observational constraint may seem 
rather weak, requiring only that (3 > 7.55 x 10 _13 m _1 , 
however it is comparable to (but larger than) existing 
constraints on similar, ad hoc, additions to GR. In par- 
ticular, constraints on additions to the Einstein-Hilbert 
action, of the form R 2 and R tlv R ,J ' v , are of the order of 
f3 R 2 > 3.2 x 10~ 9 m -1 , where (3 R 2 is the /3 parameter 
associated with the couplings of these terms [9( . Whilst 
our constraint is several orders of magnitude weaker than 
these, it will rapidly be improved as more binary pulsars 
are discovered and the observations of existing systems 
improve. This is to be contrasted with the existing con- 
straints which rely on the perihelion precession of Mer- 
cury, the accuracy of which is unlikely to improve signif- 
icantly in the future. 

As an example, white dwarf binaries reach orbital 
frequencies of the order of ~ (10 — 100)mHz towards 
the end of their merger, whilst neutron binaries reach 
~ (10— 100)Hz and should be readily observable with the 
Laser Interferometer Space Antenna (LISA) [18| ■ Whilst 
such systems would require a greater understanding of 
the strong and near field effects than that presented here, 
one can expect the constraint on j3 coming from such ob- 
jects to be of the order of j3 > (lO" 10 - 10 -6 )m -1 . 

We were able to constrain the natural length, defined 
through the /o = /(0) momentum of the cut-off function 
/ — a real parameter related to the coupling constants 
at unification — at which the noncommutative effects 
become dominant, by purely astrophysical observations. 
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